On the use of edge cracked short bend beam specimen for PMMA fracture toughness testing under mixed-mode I/II

In this paper an inclined edge cracked short beam specimen subjected to symmetric three-point bend loading was designed and examined for conducting mixed-mode I/II fracture toughness experiments. The aspect ratio (i.e. length to width ratio) and the loading span distance are considered much lower than the other conventional cracked bend beam samples. Crack tip parameters such as stress intensity factors and T-stress were computed numerically for this specimen by several finite element analyses and it was demonstrated that the specimen is able to produce full combinations of mode I and II including pure mode II. The practical capability of the short bend beam specimen was studied experimentally by conducting a set of mixed-mode fracture tests on PolymethylMethacrylate (PMMA) as a well-known model brittle material. The critical stress intensity factors, the direction of fracture kinking and the path of fracture trajectory were investigated both experimentally and theoretically using two stress and strain-based fracture criteria. The fracture toughness of tested PMMA was decreased by moving towards mode II case due to the effect of T-stress on the fracture mechanism of the short bend beam specimen.     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

Mathematical analysis of the dynamics of COVID-19 in Africa under the influence of asymptomatic cases and re-infection

Coronavirus pandemic (COVID-19) hit the world in December 2019, and only less than 5% of the 15 million cases were recorded in Africa. A major call for concern was the significant rise from 2% in May 2020 to 4.67% by the end of July 15, 2020. This drastic increase calls for quick intervention in the transmission and control strategy of COVID-19 in Africa. A mathematical model to theoretically investigate the consequence of ignoring asymptomatic cases on COVID-19 spread in Africa is proposed in this study. A qualitative analysis of the model is carried out with and without re-infection, and the reproduction number is obtained under re-infection. The results indicate that increasing case detection to detect asymptomatically infected individuals will be very effective in containing and reducing the burden of COVID-19 in Africa. In addition, the fact that it has not been confirmed whether a recovered individual can be re-infected or not, then enforcing a living condition where recovered individuals are not allowed to mix with the susceptible or exposed individuals will help in containing the spread of COVID-19.     

基于Hopf-Cole变换的Burgers方程初边值问题的Sinc-Galerkin法

利用Sinc-Galerkin法数值求解Burgers方程的初边值问题。首先,用Hopf-Cole变换将二阶非线性的Burgers方程变换为二阶线性方程,同时把第一类边界条件变为第二类边界条件。时间上的导数采用θ加权格式离散,空间导数采用Sinc-Galerkin法离散,端点处分别引入权函数处理变换后的第二类边界条件。最后,通过数值算例验证了Sinc-Galerkin法的指数收敛性,与精确解相比,本文构造的数值格式精度高,能够有效捕捉激波等物理现象。     

垂直旋转流场中单颗粒运动状态判别及受力分析

为研究单颗粒在旋转流场中的运动状态及受力情况,以毫米级球形颗粒为例,利用旋转流场颗粒运动装置,通过使用摄像机记录颗粒在流场中的运动轨迹以获取其运动参数,分析了不同转速和颗粒直径条件下颗粒的运动轨迹,拟合得到了颗粒运动状态判别公式以及颗粒运动轨迹公式,分析了颗粒在旋转流场中的受力情况。结果表明,颗粒在旋转流场平衡状态下运动状态主要分为两类,一类是未离开壁面保持静止,另一类是离开壁面保持稳定周向运动;颗粒进行周向运动的轨迹为椭圆形,并且圆心随着转速的增大靠近旋转中心,而随着粒径的增大靠近壁面;颗粒在旋转流场的运动过程中主要受到离心力和旋转科式力作用。     

On the Lipschitz equivalence of self‐affine sets

Recently Lipschitz equivalence of self‐similar sets on has been studied extensively in the literature. However for self‐affine sets the problem is more awkward and there are very few results. In this paper, we introduce a w‐Lipschitz equivalence by repacing the Euclidean norm with a pseudo‐norm w. Under the open set condition, we prove that any two totally disconnected integral self‐affine sets with a common matrix are w‐Lipschitz equivalent if and only if their digit sets have equal cardinality. The main methods used are the technique of pseudo‐norm and Gromov hyperbolic graph theory on iterated function systems.     

Collective risk models with dependence

In actuarial science, collective risk models, in which the aggregate claim amount of a portfolio is defined in terms of random sums, play a crucial role. In these models, it is common to assume that the number of claims and their amounts are independent, even if this might not always be the case. We consider collective risk models with different dependence structures. Due to the importance of such risk models in an actuarial setting, we first investigate a collective risk model with dependence involving the family of multivariate mixed Erlang distributions. Other models based on mixtures involving bivariate and multivariate copulas in a more general setting are then presented. These different structures allow to link the number of claims to each claim amount, and to quantify the aggregate claim loss. Then, we use Archimedean and hierarchical Archimedean copulas in collective risk models, to model the dependence between the claim number random variable and the claim amount random variables involved in the random sum. Such dependence structures allow us to derive a computational methodology for the assessment of the aggregate claim amount. While being very flexible, this methodology is easy to implement, and can easily fit more complicated hierarchical structures.     

The minimum forcing number of perfect matchings in the hypercube

Let $M$ be a perfect matching in a graph. A subset $S$ of $M$ is said to be a forcing set of $M$, if $M$ is the only perfect matching in the graph that contains $S$. The minimum size of a forcing set of $M$ is called the forcing number of $M$. Pachter and Kim (1998) conjectured that the forcing number of every perfect matching in the $n$-dimensional hypercube is $2^{n?2}$, for all $n\ge 2$. This was revised by Riddle (2002), who conjectured that it is at least $2^{n?2}$, and proved it for all even $n$. We show that the revised conjecture holds for all $n\ge 2$. The proof is based on simple linear algebra.     

Multistrain edge‐based compartmental model on networks

Multistrain diseases, which are infected through individual contacts, pose severe public health threat nowadays. In this paper, we build competitive and mutative two‐strain edge‐based compartmental models using probability generation function (PGF) and pair approximation (PA). Both of them are ordinary differential equations. Their basic reproduction numbers and final size formulas are explicitly derived. We show that the formula gives a unique positive final epidemic size when the reproduction number is larger than unity. We further consider competitive and mutative multistrain diseases spreading models and compute their basic reproduction numbers. We perform numerical simulations that show some dynamical properties of the competitive and mutative two‐strain models.     

Unconditionally Superclose Analysis of a New Mixed Finite Element Method for Nonlinear Parabolic Equation

This paper develops a framework to deal with the unconditional superclose analysis of nonlinear parabolic equation. Taking the finite element pair $Q_{11}/Q_{01} × Q_{10}$ as an example, a new mixed finite element method (FEM) is established and the $τ$ -independent superclose results of the original variable $u$ in $H^1$-norm and the flux variable $\mathop{q} \limits ^{\rightarrow}= −a(u)∇u$ in $L^2$-norm are deduced ($τ$ is the temporal partition parameter). A key to our analysis is an error splitting technique, with which the time-discrete and the spatial-discrete systems are constructed, respectively. For the first system, the boundedness of the temporal errors is obtained. For the second system, the spatial superclose results are presented unconditionally, while the previous literature always only obtain the convergent estimates or require certain time step conditions. Finally, some numerical results are provided to confirm the theoretical analysis, and show the efficiency of the proposed method.